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Abstract 

In this paper, we use multivariate splines to investigate the volume 
of polytopes. We first present an explicit formula for the multivariate 
truncated power, which can be considered as a dual version of the fa- 
mous Brion's formula for the volume of polytopes. We also prove that 
the integration of polynomials over polytopes can be dealt with by the 
multivariate truncated power. Moreover, we show that the volume of the 
cube slicing can be considered as the maximum value of the box spline. 
Based on this connection, we give a simple proof for Good's conjecture, 
which has been settled before by probability methods. 
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1 Introduction 

Box splines and multivariate truncated powers were first introduced in [5] and 
[8], respectively. They have a wide range of important and varied applications 
in numerical analysis and approximation theory. From the point of view of 
discrete geometry, box splines and multivariate truncated powers are closely 
related to the volume of cube slicing and the volume of polytopes, respectively. 
However, people working in the discrete geometry do not seem to be fully aware 
of the means of multivariate splines. The aim of this paper is to recast some 
problems related to the computation of volumes of polytopes and solve them 
by multivariate splines. We believe that the means of multivariate splines shed 
some light on problems concerning polytopes. 
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The main results in this paper are as follows. The exact computation of the 
volume of a polytope P is an important and difficult problem which has close 
ties to various mathematical areas. Brion's formula in continuous form (see 
[HI3)i which is well known in discrete geometry, gives an explicit formula for the 
volume of polytopes. However, the formula requires the vertex representation 
of polytopes and the generators for each vertex cone. Based on the multivariate 
exponential truncated power [19) , we give an explicit formula for the multivariate 
truncated power, which can be regarded as a dual version of Brion's formula. In 
|13j . Lasserre gave a recursive formula for computing the volume of polytopes, 
which has become a popular method today. We re-prove the formula by an 
iterative formula for the multivariate truncated power [16] , which was presented 
by Micchelli. 

Integration of continuous functions over polytopes has important applica- 
tions. For example, in most finite element integration methods, the domain of 
integration is decomposed into polytopes. Hence the integration of real functions 
over polytopes is always required. In [12) an exact formula for the integration 
of polynomials over simplices is presented. An iterative formula for computing 
the integration over polytopes is also given in [T3j. We shall show that inte- 
gration of polynomials over polytopes can be dealt with by the multivariate 
truncated power and consequently we present an explicit formula for the inte- 
gral of polynomials over polytopes. As continuous functions on a compact set 
can be uniformly approximated by polynomials, this result provides an approx- 
imate formula for integrating continuous functions over polytopes. Moreover, 
this result also shows that we can compute the integrals of polynomials over 
polytopes by calculating the volumes of polytopes. 

The volume of cube slicing is another active research topic in discrete geom- 
etry (see [52]). Suppose that Q n := [0, 1)™ is the unit cube in K n and H is an 
n — 1 dimensional hyperplane of 1" through of its center. According to [11] . 
Good conjectured that vol(H n Q n ) > 1. Hensley unexpectedly introduced a 
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probability method into the study of this conjecture and finally solved it [TT] . 
In fact, the conjecture can be reformulated as the following box spline problem: 

max-B(x|(ai,...,a n )) > , =, (1.1) 

where B(x\(a\, ■ ■ ■ ,a n )) is a univariate box spline (cf. Section 3), and are 
positive real numbers for 1 < i < n. Based on the Fourier transform of the box 
spline, we give a simple proof of (jl.ip . Hence we present a spline method for 
proving the conjecture. 

The problem of computing the volume of the intersection of Q n and an 
(n— l)-hypcrplane is also interesting and it can be traced back to Polya's thesis. 
In [18 , the authors derived a formula for the volume of such domains using 
combinatorial methods. Using box splines, we give an explicit formula for the 
volume of convex bodies obtained by intersecting Q n and a j-hyperplane, where 
j is a positive integer < n. The formula in [18) may be considered as a special 
case of ours. 

The paper is organized as follows. After recalling some definitions and nota- 
tions in Section 2, we show (Section 3) the connection between the multivariate 
truncated power and the volume of polytopes. In Section 4, we transform the 
integration of polynomials over polytopes to a problem concerning the multi- 
variate truncated power. In Section 5, we investigate the volume of cube slicing 
using box splines. Finally, Section 6 illustrates the application of the formulas 
given in this paper with some examples. 

2 Definitions and Notations 

A convex polytope P is the convex hull of a finite set of points in M. d . Throughout 
this paper, we shall omit the qualifier "convex" since we confine our discussion to 
such polytopes. Moreover, we use d-polytope to mean a d-dimensional polytope. 
When the polytope P is defined as the convex hull of a finite set of points in 
M. d , the finite set is called as a vertex representation or simply V-representation 
of P, while, if P is defined as {x € K™ | Mx = b] for some s x n matrix M and 
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s- vector b, then the pair (M, b) is called a half space representation or simply 
^-representation. For a vertex v of P, we define the vertex cone of v as the 
smallest cone with vertex v that contains P. If P C R n is a d-polytope, then let 
vol„(P) denote the d-dimensional volume of P in R™. For a rational polytope 
P, i.e., a polytope whose vertices have rational coordinates, let Rp denote the 
space that is spanned by the vertex vectors of P. The lattice points in Rp 
form an Abelian group of rank d, i.e., Rp n Z d is isomorphic to Z d . Hence 
there exists an invertible affine linear transformation T : Rp — > R d satisfying 
T(Rp (~1 Z") = Z d . The relative volume of P, denoted as vol(P), is just the 
d-dimensional volume of the image T(P) C R d . For more detailed information 
about the relative volume, the reader is referred to [5U] . 

Throughout this paper, Z + and R + denote the non-negative integer and 
non-negative real sets, respectively. Given a set D, let xd{%) = 1 if x € D, 
otherwise let Xd{%) = 0. Elements of R s can be regarded as row or column 
s-vectors according to circumstances. Let M be an s x n matrix. Then M can 
be considered as a multiset of its columns. The cone spanned by M, denoted 
by cone(Af), is the set {^2 mGM a m m\ a m > for all to}. Moreover, we set 
[[M)) := {J2m£M a ™jn\ < a m < 1, for all to g M}. Furthermore, we use #A 
to denote the cardinality of the finite set A. " • " stands for scalar product and 
|| • || for the Euclidean norm. E sxs denotes the s x s identity matrix. As a final 
piece of notation, A~ T := (A _1 ) T when A is an invertible matrix. 

3 Multivariate truncated powers and the vol- 
ume of polytopes 

Let M be an s x n real matrix with rank(M) = s. Recall that M is also 
viewed as the multiset of its column vectors. Throughout this section we always 
assume that the convex hull of M does not contain the origin. The multivariate 
truncated power T(-|M) associated with M, first introduced by Dahmen [S], is 
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the distribution given by the rule 

T(x\M)cj)(x)dx = I cj)(Mu)du, 4> G 9(R S ), (3.1) 



where £^(R S ) is the space of test functions on R s . If we define P :— {y G 
Rl\ My = x}, then (see [5]) 

T{ x\M) = - VOl " (P) (3.2) 

yj dot (MM T ) 

Note that 



vol n (P) _ v /dct(MAf T ) 



voi(P) #{[[M r ))nP} 

provided M is an integer matrix (see pQ). Hence we have 

™ M) = #{[[A^»nZ'} (3 ' 3) 

provided M is an integer matrix. In particular, if E sxs C M then T(x\M) = 
vol(P), since #{[[M T )) D Z s } = 1. So, the relative volume of polytopes P can 
be obtained by computing T(x\M). 

In the following, we shall give an explicit formula for T(x\M). We first 
introduce the multivariate exponential truncated power E c (x\M) associated with 
a complex vector c = (ci,...,c„) G C" and a matrix M. E c (x\M) is the 
distribution given by the rule (see |19j): 

E c {x\M)cj){x)dx = { exp(-c-u)cj>{Mu)du, <j){x) G (3.4) 



It is convenient for us to index the constants ci, . . . , c„ by an element m, G M, 
that is, we set c mi :— Ci for i = 1, . . . , n. For the submatrix M' = (m^ , . . . , rrii k ) 
in M we set c M ' := (c^ , . . . , c ik ) and M'/ cu> '■= {m ll /c tl , . . . , m ik /c ik ). 

We recall an explicit formula for E C (-\M). In this formula, we denote, given 
a square invertible Y C M, By ■= Y~ T cy and ay :— Y\ y eM\Y (®y • y — c y ) 1 . 

Lemma 3.1. ([15]) 

P c (x|M) = J] a y P Ci >|y), (3.5) 

rcM 

#i- = rank(y) = 3 

/or a/Z c G C™ smc/i that the denominators in ay do not vanish. 
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Using Lemma |3. II we can now give an explicit formula for T(x\M). 
Theorem 3.1. 

T(x\M) = 1 V a Y \detY\-\-e Y ■ x) n - s X con.{Y) (x) , 

[n — s)'. z — ' 

#V=iank(y)=s 

for all c € C" such that the denominators on the right-hand side do not vanish, 
where both ay and 9y are defined in Lemma \3.1[ 

Proof. For an invertible Y C M, one has 

E cy (x\Y) = , de * exp {-9 Y ■ z)Xconc(Y) (»• 
According to Lemma I3~T1 for p £ M \ 0, we have 

E pc {x\M) = p- n+s a Y E pCY (x\Y). 

Y<CM 
#y=rank(y)= S 

Then the Taylor expansion of E pc (x\M) about in the variable p is 

oo 

E pc (x\M)=J2p l ~ n+S PlW> ( 3 ' 6 ) 
1=0 

where 

= 7? X! ay det F| -1 (-0y • a;) Z Xcono(y)(a;)- 

y c m 

#Y=i.nk(Y) = s 

The definition of T(a;|M) implies that it is the constant term in (|3.6[) . Hence, 
we have T(x|M) = p n - s (x). The theorem follows. □ 

Brion's formula, which is obtained by Brion's Theorem, is useful for com- 
puting the relative volumes of polytopes. We state it here. 

Theorem 3.2. ([7]) Suppose thatV is a simple rational convex d- poly tope. For 
a vertex cone tC v of V , fix a set of generators wi(v), W2(v), . . . , Wd(v) 6 Z d . 
Then 

(-if (?;■ c) d | det( Wl (?;),..., 

V a vertex of V Uk=l ( W k ( v ) ' c ) 



/or c £ C d suc/i that the denominators on the right-hand side do not vanish. 
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Remark 3.1. Brion's formula requires the V -representation and the generators 
for each vertex cone, while the formula presented in Theorem \3.1\ requires the 
Ji-representation. Hence, the formula in Theorem \3.1\ can he considered as a 
dual version of Brion's. 

We next turn to another formula for computing the relative volume of poly- 
topes. We first introduce an iterative formula for calculating the multivariate 
truncated power. 

Theorem 3.3. (|16j) Let M be an s x n matrix with columns mi, . . . ,m n € 
R s \{0} such that the origin is not contained in conv(Af). Suppose that n > s+1. 
For any Ai, . . . , A„ £ K and x — Sj=i ^j m j> we have 

1 " 

T(x\M) = >>jT(x\M \ rrij), (3.7) 

n ~ S i=i 

In [13j , Lasserre found a formula which expressed the volume of polytopes as 
the combination of the volumes of the faces. We describe it as follows. Consider 
the convex polytope defined by 

D(b) := {x e R"| Ax < b}, 

where A is an s x d matrix and b an s-vector. The ith face of D{b) is defined by 

Di(b) := {x e R d \ a t ■ x = b u Ax < b}, 

where ai is the ith row of A. We set V(d, A, b) = vo\d(D(b)) and set Vi(d — 
l,A,b) = vo\d(Di(b)). Now we can describe Lasserres formula and present a 
proof by (|3~7)l . 

Theorem 3.4. ([13]) IfV(d, A, b) is differ entiable at b, then 

V{d,A,b) = -j2^rMd-l,A,b). (3.8) 

Proof. Without loss of generality, we can suppose that all points in D(b) are 
non-negative, i.e., D(b) :— {x £ Ri_\ Ax < b}. We first consider the case where 
each entry in A is an integer. By (13. 2|) and (|3.3|) . when A is an integer matrix, 

T(b\M) = vol(P) = V(d,A,b), 
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where P := {x G M+ +s | Mx = b} and M := (A,E sXs ). Let be the s-vector 
with 1 at the ith position and for j ^ i. Using (|3.7|) . we obtain 

T(b\M) = \Y,bi T(b\M \ e { ). (3.9) 

i=l 

Note that the (d— l)-polytope Di(b) lies in a hyperplane {x G M. d \ a^x = bi} and 
that the (n— l)-dimensional volume of the unit parallelogram in the hyperplane 
^ s gcdfJ^'^ a ) CP)' wnere a ij i s the jth entry in the vector a^. So, one has 



vol d (A(6)) H| 



vol(A(&)) gcd(a il ,...,a irf )' 
By (J331) and (I^TO]) . we have 

vol(A(6)) vol(A(6)) 



(3.10) 



T(&|M\e;) 



#{[[M\eO)nZ s } gcd(a i i,... ! a«) 
void (A W) ^(d-l.Afe) 
l|Oi|| lloiH 



Substituting T(b\M) = V(d,A,b) and T(6|M \ e 4 ) = ^gjp^ into (l3~9l) . 
we get (|3 . 8[) . By taking limit, ()3 .8[) holds for any matrix A. □ 

4 Integration of polynomials over polytopes 

In this section, we consider the problem of integrating of polynomials over 
polytopes. Since each polynomial can be written as the sum of monomials, 
we only consider the monomial case. For every k = (ki, . . . ,k n ) G and 
M = (mi, . . . , m n ) we set 

M k := (mi, . . . ,TOi,m 2 , . . . ,m 2 , . . . ,m n , . . . ,m n ). 



fel+l fe 2 +l fen+l 

The following theorem shows that the integration of monomials can be handled 
by the multivariate truncated power. 

Theorem 4.1. Suppose k = (k\, . . . ,k n ) G Z" and f(u) = YYj=i M j J • 
P := {u G R!^| Mu = a;}. T/ien 

T(a;|Af k ) = 1 / f(u)du, 

k! • \J dct(MM T ) 

where k! := fci! • • • fc ra !. 
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Proof. We set 



T k (x\M) := 1 j f(u)du 



and consider its Laplace transform, i.e., 

Tk(u\M) := I exp (-lo ■ x)T k (x\M) dx. 



For every u in P := {u G M™ Mm = x}, we can write u in a unique way as 
u = z + y where z £ ker(M) and y_Lker(M). Note that x = My and hence 



y/ ' det(MM T )dy = dx (see [5]). Then we have 
/ exp (— u> ■ x)Tk(x\M)dx 



1 



exp (—lo ■ x) I f{u)dudx 

y>det(MM T ) Jr* ./{mGR™ Mu=x} 

exp (-lo ■ My) / x K » ( z + v)fi z + y)dzdy 

yJ_kerM J zEkerM + 

/ x K « i z + y)f( z + y) cx p (~ u ■ M ( z + y)) dzd v 

f - 1 

/ f{u) exp (-to ■ Mu)du = k! TT -7—3-. 

, i^l v - // 



Also, note that 



™ 1 

T(^|Af k ) = ff 7 

( u . m j ) k 3+ 1 



3- 

Hence the theorem follows from the inverse theorem for Laplace transform. □ 

5 Multivariate box splines and the volume of 
cube slicing 

The multivariate box spline B(-\M) associated with M is the distribution given 
by the rule (see [H[S]) 



B(x\M)(j)(x)dx = / <j)(Mu)du, e @(R S ). (5.1) 

[0,1)" 



According to [6] , one has 



b(x\m) = X£U££jW (5 . 2) 

1 ' ' ^det(MM T ) 1 ' 
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where P := {y E R™| My = x}. The formula (|5.2[) shows the connection 
between the box spline and the volume of cube slicing. Based on this connection, 
we can study some interesting problems concerning the unit cube. 

Recall that Q n is the unit cube in R". We define H is an n — 1 dimensional 
hyperplane of W 1 through the center of Q n , i.e., 

m 

H := {y e R™ am + ■ ■ ■ + a m y m = ^a,/2}, 

i=l 

where 1 < m < n , a, is a real number for 1 < i < m. Based on the symmetry 
of Q n , we can suppose > for 1 < i < m. We set A :— (ai, . . . , a m ). Then 
from (15.21) we have 



vol n {H R Q n ) = y det(AA T )B ((oi + • • • + a m )/2|A) . 

By the symmetry of the box spline, -B(x|j4) achieves its maximum value at 
(di H h a m )/2. So, as stated before, Good's conjecture is equivalent to 

maxB^IA) > * (5.3) 
V Z^i=i a i 

where 1 < m < n. We next present a spline method for proving 
Theorem 5.1. 

maxB(x|(oi, . . . , a m )) > 



where a; are positive real numbers. The equality holds if and only if m = 1. 
Proo/. Set 

m 

C(ai|A) := B(x + ^a l /2\A). 

i=l 

The Fourier transform of C(x\A) is 

?— l ' 

We can see that C(0|^4) = 1. According to the definition of Fourier transform, 



we conclude that 



° 1 - V m a 2 

t 2 C(t\A)dt = -~C"(0\A) = (5.4) 
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Put 

S(t) := f C(x\A)dx. 
Jo 

By 

(tmaxC(x\A)) > S(t) 2 



we have 



(maxC(x\A)j j t 2 C(t\A)dt > j S(t) 2 C (t\A) dt (5.5) 

1 [°°dS(t)Z 1 



3 7 dt ~ dt 24' 

where the last equality follows from 

' C(0\A) 1 

C(x\A)dx = -±±-L = -. 

We combine f|5.4[) and (|5.5[) to obtain 

max B(a;|A) = ma,xC(x\A) > 



vE,=i«, 

From (|5.5|) , we see that the equality holds if and only if m = 1 . □ 



In the following theorem, we present an explicit formula for the volume of 
j-slice of Q n . The formula given in [TB] is the (j = n— l)-case in the following 
theorem. 

Theorem 5.2. Suppose that M is an (n — j) x n matrix with rank(M) = n — j 
and let P := {y G R£| My = x}. Denote S = {0, 1}™ and \e\ = £" =1 £ *- Then 
we have 



. / r~t /~\ \ \/ det(MM T ) 

voi„ (P n Q n ) = 2 K - '- 

r- 

J2 <xy\ det(y)!- 1 J2(-^ s1 (-0y • (x - Me)y X con C (Y)(x - Me), 

YCM £ gH 
#y = rank(y) = s 

for all c such that the denominators on the right-hand side do not vanish, where 
both ay and 9y are defined in Lemma \S.l\ 
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Proof. By (|5.2p , one has 

vol n (P n Q n ) = yJdet(MM T )B(x\M). 

Now we present an explicit formula for B(x\M). Recall that B(x\M) = V mT(x\M) 
(see 0), where V M := U7=i v ™* and V m4 T(ai|M) = T(x|M) - T(ai - m !; |M). 
By the formula and the property of the difference, the following formula can be 
obtained (see [T7]): 

B(-\M) = ^(-l)l £ lT(- - Ms\M). (5.6) 
The theorem is proved by putting (|3.6|) into (15. 6j) . □ 

6 Examples 

Example 6.1. 5e£ 

£>(«) := {y e j/ x + y 2 < z; -2y x + 2y 2 < z; 2y x - y 2 < z}. 

The volume of D(z) has been calculated in \15]j using Cauchy's Residue theorem. 
Here we can obtain it directly. Based on I13.3\) . we have vol(£>(z)) — T(z\M) 
where 

/ 1 1 1 0\ 

M = -2 2 1 , z = (z, z, z) T . 
\2 -10 1/ 

We use mi to denote the ith column in M . A brief calculation shows that the 
cones spanned by the square matrices (mi, m 2 , m^), (mi, m 2 , m§), (mi, 7773, m^), 
(7772,7773,7715) and (7713,7714,7715) contain z. We select c — (1,1, 1,1,1/2) in The- 
orem HOI and obtain that T(z,\M) — j^z 2 which agrees with the result presented 
in 115$. 

Example 6.2. Set fid := {y G WL\ Yli=i Hi — !}■ We consider the problem of 
integrating of monomials over fid, i-e., 

Jd ■= I y k i---y k d d dy x ---dy d . 
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The value of Jd is also calculated in fl(A \21]j . Based on Theorem \4-l\ we can 
compute it easily. We set := (1, . . . , 1) € 7L d . Using Theorem \4-l\ we have 



Jd = fci!- •• k d l 



T(l\e. 



) 



.=1 



It is well known that T(x\ed) 




- (d-iy. 



Hence, we have 



J d = kx\---k d \ r(l|e^ 




hi- ■ -k d \ 
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